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We determine the numerical values of scalar multi-loop two-vertex Feynman diagrams, the gen- 
eralized sunset diagrams, by integrating all but the longitudinal momenta analytically. For the 
longitudinal momenta we introduce one collective coordinate, which allows us to determine the nu- 
merical value of the diagram efficiently and to an arbitrary accuracy. The imaginary part and the 
threshold behavior of the diagram are also handled within this framework. 
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I. INTRODUCTION 

Multi-loop Feynman diagrams pose a serious challenge, especially in the massive case where multiple scales arise. 
Apart from the successes of high-order expansions in a in QED, which provide a stringent test of quantum field 
theory, Feynman diagrams are also the way in which we understand the most of field theory, since it allows us to 
quantify it. They form the basis for our understanding of many phenomena, such as asymptotic freedom and gauge 
invariance. Therefore, we should strive to refine our handling of Feynman diagrams and to extend its context and its 
applications. 

In the late 40's, the move from time-ordered perturbation theory |]l[ towards covariant perturbation theory brought 
about a revolution in field theory. |2| It created a handle on the calculations and people were able the control the 
divergences. However, despite the successes for scattering experiments, the method failed to deliver for bound- 
states calculations. Therefore, there has been a constant move "backwards" to quasi-potential, and time-ordered, 
formulations, as we can formulate a bound state only in a single time frame, not with the relative time as it exists 
in the covariant formulation. However, often it is hard to follow the route back from covariant to time-ordered 
perturbation theory, and then to apply it to an intrinsic non-perturbative problem, such as a bound-state problem. 
One distinct method for describing bound-states in field theory is discrete light-cone quantization, where as the time 
direction a light-like direction is chosen. This has certain advantages, extensively discussed in the literature. j!| 
However, also here renormalization forms a serious problem, often dealt with rather callously. Although the renor- 
malization for perturbative expansions is reasonably under control, Q,|5| its extension to non-perturbative calculations 
is far from desired. 

In this paper I will tackle a class of simple n-loop Feynman diagrams, determine the finite part and show that the 
light- front approach is particularly useful for that. I will also discuss how to look upon the large set of counterterms 
in this highly divergent case. These diagrams have been studied extensively in the recent years. |(| However, the 
simplicity of this approach in Minkowski space is striking. There are no special functions needed, and eventually it 
depends on the introduction of the collective coordinate (3, analogous to the radial coordinate r 2 in Euclidean space. In 
this case, specifically, the coordinate interpolates smoothly from the threshold value at the center of the kinematical 
to the edges of the kinematical domain. Eventually, such simple collective coordinates for many-particle systems 
might extend the applicability of Hamiltonian light-front field theory, as, generally, the problem of the Hamiltonian 
approach is the control on number of the variables. 



II. THE SUNSET DIAGRAM 



We consider the n-loop Feynman diagram, which consists of n + 1 lines between two vertices: 
1 f d 4 k 1 ---d i k n 



(2™)" J (kl - m\){kl -ml)--- {kl - ml)((p -k 1 -k 2 k n ) 



(2.1) 



*e-mail: ligterin@ect.it 



1 



This diagram the generalized sunset diagram. Lines can be added or removed (see Fig. ||). The Feynman diagram 
is covariant and therefore it is only a function of p 2 and the masses mi, • • • , m n _|_i. We solve it in the frame where 
p± = (p 1 ,p 2 ) — 0. We introduce light-front coordinates: fc ± = -^(fc°±fc 3 ), and Xi = kf /p + . The transverse momenta 

k 1 and k 2 are unaltered. After the residue integration over the light-front energies k~ , we find: |7) 



where 



a= ^ + --- + ^+ ^ J;" . (2.4) 

X\ X n 1 %i 

The domain A is given by Xi > and ^2 x i < 1- This integral is the corresponding light-front diagram, equivalent to 
the Feynman diagram. If we translate the transverse momenta successively starting from k± n : 



i±<~ \Tk± 3 r - + — =5 — = , (2.5) 



l n = I cfcx • • • d,„ / d^ • • • d^ ^.,.^^^^,^,^ , (2.2) 



/ r 1 = ^--- + ^ + T ^L_ , (2.3) 
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we find that the a reduces to a pure quadratic form in 

(2.6) 



The integral is divergent, and requires counterterms Co, cip 2 , • • •, c n _i(p )™ . We do so by subtracting the (n— l)-th 
order Taylor expansion in p~ around p 2 — , implemented through the the multiplication of the integral with the 
proper moments: Q 



where 



which will lead to 



J n = I -> Z£* = / 7J n , (2.7) 

J A J A 



v 2 

J=l^— ■ (2-8) 



rreg _ 



0- 1 + a 
dxi ■ ■ ■ dx n d 2 l±i ■ ■ ■ d 2 l± n (p 2 ) n 



I a 2" +1 £i ■ ■ ■ x n {l - E«i)(2p + P- - P- 1 - a)^- 1 + a) n 
If we scale the transverse momenta accordingly and perform the angular integrations we obtain the integral: 



(2.9) 



00 r f^2\n 



where 



1 ,hi -i 1 "'•■•^PFri^rw ' (2 ' 10) 



Zi = ll.(l + --1 ] . (2.11) 



3=1 ^ 



In order to integrate over Zi we express the integrand as a series in p 2 . For each separate term we can integrate over 
all Zi, which yields: 

F°g = --?—/ d"^^ 2 )"- 1 Y - ( P 2 f3) l+1 . (2.12) 

2 n + 1 J A v ' (n + i) ■ ■ ■ (i + 1) K ' v y 
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We can write the series as an analytical function: 



X 2^(n - l)!/3n-i (1 ~ T "- l} ( X ~ ln " P ^ ~ ^ ' (2 ' 13) 

where l~ n -i stands for the (n— l)-th order Taylor expansion around p 2 = 0. The Taylor expansion yields the following 
polynomial: 

n—l I k (—T^j ! \ 



The imaginary part follows directly from the natural logarithm of Eq. (2.13) 



lim ln(— x — ie) = ln |x| — ii{6{x) , (2-15) 



which lead to a finite amplitude, unaffected by the renormalization procedure which reminiscence appears in Eq. ( 2.13| ) 
in the form of the subtracted Taylor expansion. 



III. LONGITUDINAL INTEGRATION 



After the integration over the k i 's and the fc_u's, we are left with an integration of the lo ngitu dinal momentum 



fractions Xi, which cannot be performed analytically. However, note that the integrand of Eq. (2.13) depends only on 
one particular combination of the longitudinal momenta, namely (3. This j3 is a smooth function of the longitudinal 
momenta Xi, and ranges between: 

0<P<\J2m?J =b , (3.1) 

where for (3 — b the longitudinal momentum fractions Xi equal rriiVb. Therefore the n-dimensional integral over A 
reduces to the determination of the integration measure /i for the integration over /3: 



d n xf{fi) = / KPWf(P) ■ (3.2) 

A JO 

The volume of the domain A is T(n + 1) . Once this measure is determined, it can be used for all values of p 2 . The 
measure fi can be determined via several means, for example, with Monte Carlo integration. The threshold behavior 
of the diagram is dominated by the the values of f3 close to b. For this purpose we can make an analytical expansion 
of around b. We find that: 



n+l 
.i=l 



(6-/3)-i- • (3.3) 



where ft n is the surface area of the unit sphere in n dimensions. Note that as some of the masses tend to zero, the 
exponent in the measure will be larger than (n — 2)/2. The addition of a zero mass particle, m.j = 0, leads to a 
flat direction in (3 with respect to the longitudinal momentum fraction Xi at the threshold. Therefore the harmonic 
approximation breaks down. However, [3 and the measure fi are well-defined as long as at least one particle is massive. 

A. The integration measure 

Apart from series expansion and Monte-Carlo integration mentioned above, we can determine the measure iteratively. 
Given the integration measure /U n (/3 n ) for n momenta, the integration measure /i n _|_i(/3 n _|_i) for n+l variables can be 
expressed as 



bn+i rb n rl 

fhi+i{(3 n +i)d{3 n+1 = I / iJ. n (Pn)df3 n y n dy , (3.4) 
Jo Jo 
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where 
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1 - V PnV 



(3.5) 



with m the mass of the added particle, with longitudinal momentum fraction 1 — y. The other longitudinal momenta 
are scaled by a factor y such that the total longitudinal momentum remains 1. 



In this paper we described a method to find the finite part of any diagram of the type of Fig. |lj. The divergences are 
removed using the Taylor expansion in the external momentum. The counterterms, Co, Cip 2 till Cn—i{p ) n , for the 
sunset diagram can be expressed as divergent integrals depending on the masses mi, • • • , m n +i: 



which defines their relation with other renormalization schemes. The other renormalization schemes will find that 
the counterterms in Eq. d4.1| ) equal an infinite constant, to be removed, and a finite part, a function of the masses, 
which is the finite renormalization. The use of Taylor expansion became in disfavor, because of two complications. 
Firstly, in the case of multiple external momenta, it is not clear which combination of external momenta should serve 
as variable in the Taylor expansion; different choices will lead to different results, and do not automatically guarantee 
locality. Secondly, in the case of gauge theories, an extremely consistent scheme, which treats a whole class of integrals 
in the same way, is required such that the gauge invariance is preserved. Dimensional regularization has for a long 
time been the only scheme satisfying this consistency, which preserved algebraic relations existing among different 
integrands of Feynman integrals. For example, the fermion-loop correction to the gauge propagator must be 
transverse, therefore the two parts to namely g^H s and p^p^ILt must be handled in the same way such that 
p 2 Iit — — n s , which is difficult problem for an arbitrary regularization scheme, since the two terms have different 
degrees of divergence. 

However, for an Hamiltonian approach, such as light-front field theory, the renormalization at the level of the integrand 
is required, if one wants to carry the renormalization procedure over from the covariant renormalization. Note that 
a straightforward cut-off procedure breaks covariance, since it cannot be applied to the energy part of the covariant 
integration. The integration of the energies and the regularization should be interchangeable, such that locality is 
guaranteed. So, although it leads to further complications which requires careful analysis, the Taylor expansion is the 
way forward for the Hamiltonian approach. [f| 

The natural choice of renormalization for a Hamiltonian approach, is to make the self-energy contributions vanish as all 
the particles are on-shell. However, this is not consistent with the covariant, local and therefore true, renormalization. 
If the sum energy is the sum of the on-shell energies, it does not mean that the energy is shared out evenly; a large 
part of the amplitude might arise from the case that both particles are off-shell in different directions. Therefore it is 
essential to treat the subtractions as pure constants. 

Even more, although we can generate finite terms in light-front perturbation theory, for a proper light-front approach 
we should take the procedure one step further and determine the corresponding finite wave function. However, this 
is far beyond the scope of this paper. 



Central to this approach is the actual shape of the measure For particles of equal mass, we find that the 

measure is most spread over the whole range of (3. As the masses start to deviate the measure peaks more and more 
at low values of f3. However, the measure stays finite, even for massless particles. In Fig. |^ we show two scaled, 
normalized set of measures, one for equal mass particles, and one for particles with increasing masses. Note that the 
increasing masses peak more at lower values of /3, due to the leading contributions from the heavy particles carrying 
large momentum fractions. In Fig. |^ we compare the measures for two massive particles and a number of particles 
with equal, but small, or vanishing, masses. 

For the inspection of the amplitudes, I fitted the measures with a five parameter function, which fits the measure 
with an accuracy within a few percent. The accuracy for the massless case is higher than for the massive case, in the 
former case it is below a percent. The function depend on the parameters 71, 72, • • • , 75: 



IV. RENORMALIZATION 




(4.1) 



V. RESULTS 
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M(£)= 7 i(i-£r 2 +73(i-£r 5 £ 74 , (5.1) 

where = 0/b and \x = T(n+l)/j,/b, such that the axis and the measure are normalized to unity. The fitted parameters 
for the two extreme cases, one case with all masses equal.and one case with the first two masses 1.0, and all the other 
masses zero are given respectively in Table Q and Table y. 



VI. CONCLUSIONS 



I have derived a straightforward, and largely analytical, way to determine the finite part of the sunset diagram. 
Both the threshold behavior and the full amplitude can be determined accurately. I removed the divergent parts 
by subtracting the corresponding Taylor expansion. The mass dependence of the diagram, for the scaled external 
momentum bp 2 , is only weak over the whole range of masses, and this dependence appears solely in the integration 
measure 

A careful analysis of the Feynman parameterization || could also yield a similar variable 0. However, it requires one 
to work in Euclidean space, where the imaginary part does not come for free. Also the poles in dimension space, 
which are the different subdivergencies of the integral, are transferred to singularities in the parameter space, which 
renders dimensional regularization invalid. Therefore the removal of the lower-order Taylor expansion is essential. 
It seems possible to extend the light-front approach to more complicated diagrams, which could contain two and 
more light-front intermediate states. This requires the introduction of more variables. The determination of the 
measures stays essentially the same. Eventually, the light-front approach might be more convenient for the calculation 
of multi-loop diagrams, as it sees such a diagram as a transition via collection of successive intermediate states, which 
can all be handled separately, and do not grow as wildly as the number of the subgraphs of a complicated covariant 
multi-loop Feynman diagram. 
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FIG. 1. The generalized sunset diagram 




FIG. 2. Two sets of integration measures, all normalized to unity, left with one to six loops, where all the masses are 
equal, and right, with one to six loops including heavier and heavier particles: {1.0, 1.0}, {1.0, 1.0, 2.0}, {1.0, 1.0, 2.0, 3.0}, ■ ■ ■, 
{1.0, 1.0, 2.0, ■ ■ ■ , 6.0}. The variable /3 is scaled such that it runs from zero to one. 
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m=1 .0,1 .0,x,...,x 




FIG. 3. Two sets of integration measures, all normalized to unity, left with one to six loops, where the first two masses are 
1.0, and the other masses are small (0.01) for the solid line, or zero, for the dashed line. Note the rapid decline at /3 = 1. 



TABLE I. The fitted parameters for the equal mass case. 



m 




7i 


72 


73 


74 


75 


1.0,1.0,1.0 
(1.0) 4 
(1.0) 5 
(1.0) 6 
(1.0) 7 




.661422 
.746938 
.817047 
.874748 
.921155 


1.01707 
.476909 
1.08804 
1.55609 
1.94665 


1.23111 
1.66776 
3.45241 
5.90012 
9.09230 


.815085 
.814929 
.950327 
1.06571 
1.16310 


.0064803 
.526326 
.982281 
1.44554 
1.88020 






TABLE II. The fitted parameters 


for the massless case. 






m 


7i 


72 


73 


74 


75 


1.0,1.0,0.0 
(1.0) 2 , (0.0) 2 
(1.0) 2 , (0.0) 3 
(1.0) 2 ,(0.0) 4 
(1.0) 2 , (0.0) 5 


.988223 
1.48154 
1.97032 
2.45665 
2.94061 


.538897 
1.33267 
2.19627 
3.05772 
3.92001 


1.07238 
2.60100 
4.63872 
7.09066 
9.91789 


.706459 
.706458 
.702103 
.698004 
.694266 


.490497 
1.66423 
2.75778 
3.84503 
4.92679 
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